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°^ : Abstract 



o 

We investigate the contribution of zero-point motion, arising from fluctu- 
ations in kaon modes, to the ground state properties of neutron star matter 
5_i ■ containing a Bose condensate of kaons. The zero-point energy is derived via 

the thermodynamic partition function, by integrating out fluctuations for an 
arbitrary value of the condensate field. It is shown that the vacuum coun- 
terterms of the chiral Lagrangian ensure the cancellation of divergences de- 
pendent on fi, the charge chemical potential, which may be regarded as an 
external vector potential. The total grand potential, consisting of the tree- 
level potential, the zero-point contribution, and the counterterm potential, is 
extremized to yield a locally charge neutral, beta-equilibrated and minimum 
energy ground state. In some regions of parameter space we encounter the 
well-known problem of a complex effective potential. Where the potential is 
real and solutions can be obtained, the contributions from fluctuations are 
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found to be small in comparison with tree-level contributions. 



Typeset using REVT^X 



I. INTRODUCTION 



In recent years the idea that the ground state of nuclear matter at high density contains 
a Bose condensate of kaons has received considerable attention. The existence of this novel 
state of matter is suggested by the attractive s-wave kaon-nucleon interactions found in 
low energy effective chiral Lagrangians, which, by construction, preserve the symmetries of 
QCD. This is supported experimentally by kaonic atom data which indicate Q a strongly 
attractive K~ optical potential. Attractive kaon-nucleon interactions imply that in dense 
matter, such as that in the interior of neutron stars, the effective kaon mass is lower than 
in free space. It could thus be energetically favorable for leptons in the matter to convert 
to kaons via strangeness changing interactions, in which case the lowest energy kaon mode 
becomes macroscopically populated, forming a condensate. 

Since the original suggestion by Kaplan and Nelson f2|, studies of kaon condensation 
and kaon propagation in the nuclear medium can broadly be divided into three categories: 
studies of the basic interactions, consequences for heavy-ion collisions and astrophysical 
implications. 

A number of studies have been aimed at improving the description of kaon-nucleon 
interactions over that of the original model. Brown, Lee, Rho and Thorsson || included all 
terms to next-to-leading order (0(q 2 )) in chiral perturbation theory, and fitted the expansion 
parameters to KN scattering lengths. Lee, Brown, Min and Rho ||] included all terms to 
next-to-next-to- leading order (0(q 3 )) in chiral perturbation theory and also fitted kaonic 
atom data. Lee et al. have also included the A(1405) (which dominates the low-energy 
K~p interaction), four-baryon interactions, and Pauli blocking. Other noteworthy studies 
are those of Ref . || , concerning the effect of kaon-nucleon correlations, and the work of the 
Kyoto group ( see e.g. Ref. [|] ). Waas, Kaiser, and Weise have recently investigated 
kaon modes in matter using a multichannel K-matrix analysis. For our purposes, it is 
worth mentioning that the majority of these more refined investigations have shown that 
the basic mechanism for kaon condensation, namely the existence of an attractive interaction, 
is robust. 

As regards heavy ion collisions, strangeness is conserved on the relevant timescales and 
a kaon condensate is not likely to form. However, modification of the kaon modes follow- 
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ing from the chiral Lagrangian could have consequences in kaon flow ||, for subthreshold 
production of kaons 0, or for the 0-meson peak in dilepton spectra [[II]]. 

Another major direction that investigations have taken is to study the astrophysical 
consequences of the formation of a kaon condensate (the effects appear to be reduced if 
hyperons are also present |TT|, a possibility that we do not consider here). The presence of 
such a condensate considerably enhances the proton concentration which leads to a softer 
equation of state (EOS) and a lower maximum mass for the cold neutron star. The actual 
value of this maximum mass is dependent on the nuclear EOS, but for the non-relativistic 
EOS employed here it is quite close p|, [T2| . |l3| to the lower limit of 1.44M Q set by PSR 
1913+16 ||14|| . An interesting effect occurs when the neutron star has just formed since 
electron neutrinos are trapped for about 10-15 s (also the entropy/baryon is ~ 1 — 2, but 
this plays a lesser role). The changes in the stellar composition induced by trapping increase 
the threshold for kaon condensation. This stiffens the EOS and allows a maximum mass 
which is ~ 0.1 — O.2M larger than in the absence of trapped neutrinos. Thus there is 
the possibility of metastable neutron stars which become black holes after about 10-15 s 



when the trapped neutrinos have left [|13|jl5| pT8|1 . This could account for the fact that fewer 
pulsars are observed in supernova remnants than expected from conventional scenarios [ITS 



As regards dynamical simulations [15|20| of the delayed collapse of hot neutron stars into 
black holes, Baumgarte, Shapiro and Teukolsky |21] have recently examined the effect of kaon 
condensates. Here the cold EOS is supplemented with simple assumptions about corrections 
for finite temperature. 

One-loop calculations to date are valid only up to the density at which the condensate 
forms where the condensate amplitude is vanishingly small. In order to go beyond threshold 
we analyze the fluctuations around the kaon condensate to one loop order and include the 
full non-linearities of the Lagrangian. This is of general theoretical interest, as well as in 
the particular context of neutron stars. While we shall restrict our numerical analysis here 
to the zero-point contribution at zero temperature, the formalism also provides a consistent 
treatment of the thermal contributions due to kaons. These would be needed in studies of 
newborn neutron stars. 

In the work which follows, we shall make the following approximations. Nucleons are 
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treated non-relativistically, whereas a more complete treatment would involve a relativistic 
treatment to one loop order. In the meson sector, we restrict the discussion to kaonic degrees 
of freedom. We thus omit from the discussion fluctuations in the pion and eta sectors. 
Physically one might argue that ignoring non-kaonic pseudoscalar degrees of freedom is 
justified since one expects that the corresponding modes are less modified in nuclear matter 
than those of kaons. The methods developed here could be used to include the full octet 
of pseudoscalar fluctuations, but this would further complicate the present analysis. Lastly, 
we have sought to treat the chiral Lagrangian at the most simple level starting from the 
Lagrangian of Kaplan and Nelson ||, although one would wish to include the additional 
terms of Refs. at a later stage. 

The paper is organized as follows. In Sec. II, we present the basic formalism, including 
the evaluation of the kaon partition function, a discussion of divergences in the zero-point 
energy, the nucleon Hamiltonian, and the equilibrium conditions. In Sec. Ill we compare 
numerical results in neutron star matter with and without zero-point fluctuations. We dis- 
cuss the critical density for kaon condensation, the composition of the star and the EOS. 
Following the Conclusions in Sec. IV are three Appendices. Appendix A concerns the can- 
cellation of divergences dependent on the chemical potential in a toy model, while Appendix 
B discusses the case of the actual chiral Lagrangian. In Appendix C we derive the critical 
density from the kaon propagator to one-loop order, and show that it is the same as that 
obtained by expanding the zero-point energy in terms of the kaon fields. 

II. FORMALISM 

Our principal interest is the evaluation of the impact of the kaon zero-point energy 
contribution on the structure of neutron stars. In Subsec. A we obtain the partition function 
for the chiral Lagrangian of Kaplan and Nelson || to one loop order and in Subsec. B 
we discuss the elimination of the divergences. For the nucleons alone we use a simple 
non-relativistic Hamiltonian which should be sufficient for present purposes, as outlined in 
Subsec. C. The stellar equilibrium conditions that need to be satisfied are given in Subsec. 
D. 
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A. Kaon Partition Function 



The SU(3) x £77(3) chiral Lagrangian involving the pseudoscalar meson octet and baryon 
octet is MJZ 



(i) 

ct 



C K = £« + + £ff + df + 



£« = iTrBv ■ dB + zTrS^K, B] + 2/JTr77£ i/ {A y , 5} + 2FTrBS u [A l/ , B] + C 
C (2) = ^fr Trd ^ ud u U+ CTrm q (U + U ] - 2) + ai Tr B(£m q £ +h.c.)B 

+a 2 TrBB(£m q tl + h.c.) + a 3 {Tr/7/7} {Tr(m 9 [7 + h.c.)} + £ (c>d>/) + £ 



(2) 
ct 



The bracketed superscripts denote the order in the low momentum expansion defining the 
effective field theory. The various one-loop counterterms, with subscripts ct, are discussed 
in Appendices A and B; here in the main text we will simply put them all together in 
a single £ ct . The additional terms denoted by C^ c ' d ^> are given in Refs. [|3],f4| and those 
with coefficients d can contribute to s-wave kaon-nucleon interactions. In the interests of 
simplicity we set C^ c ' d '^ = and so use as a starting point the Lagrangian of Kaplan and 
Nelson [0. The Lagrangian ([I]) is derived in the limit of heavy baryons, where v v = (1, 0) 
is the velocity four- vector of the baryons, and £^ is the spin operator (v ■ S — 0,£ 2 = — |). 

In Eq. ([[]) , U = £ 2 is the non-linear field involving the pseudoscalar meson octet, from 
which we retain only the contributions- 

/ K+\ 



U = exp (^Y^ 



M 




VTT- / 



(2) 



and from the baryon octet we retain just the nucleon contributions 

/0 p\ 

B= n . (3) 
\0 0/ 

Consequently the p-wave interactions in £W with coefficients D and F do not enter (A u = 



ct£}). In Eq. flU), the quark mass matrix 



diag(0, 0, m s ); i.e., only the mass of 



the strange quark is taken to be non-zero. For the mesonic vector current, V u , only the time 
component survives in an infinite system with Vq = , <9q£]. Also, the pion decay constant 
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/ = 93 MeV, and C, a\,a 2 , and 03 are constants. After some algebra, the relevant part of 
Ck can be written 

C K = \(l + 7 )d v K + drK- + \{K + d v K-f + {K-d u K + ) 



+11 



ibK + d Q K~ - (m% + c)K + K 



C ct ■ (4) 



Here, the kaon mass is given by m? K = 2Cm s /f 2 , and we have employed the following 
definitions: 

~ 2 2 /sin6\ 2 /sin±9\ 2 



[1 + x) — ; c = (aix + a 2 + 2a 3 )— — . (5) 



p ' ' V e J ' V ie 

ft Tfl 

I— ; c = (aix + a 2 + 2a 3 )-^ 

Expectation values of nucleon operators are evaluated in the mean-field approximation. The 
total nucleon number density is n = n p + n n and the proton fraction is x = n p /n. Note 
that since we shall use a non-relativistic model for the nucleons we have not distinguished 
here between the scalar density and the number density in the expression for c. The nucleon 
kinetic energy in will be considered in Subsection C and so is not included in Eq. (f|). 

We have not included in Eq. (f|) terms which simply give a constant shift to the baryon 
masses; they indicate that aim s = —67 MeV and a2m s = 134 MeV, using the hyperon- 
nucleon mass differences. The remaining constant a^m s is not accurately known, and we 
shall use values in the range —134 to —310 MeV corresponding to to 20% strangeness 
content for the proton. The corresponding range for the kaon-nucleon sigma term, 



^KN 



-|(ai + 2a 2 + 4a 3 )m s , (6) 



is 167-520 MeV. Some guidance is provided by recent lattice gauge simulations f23| , which 
find that the strange quark condensate in the nucleon is large, i.e., (A^|ss|A^) = 1.16 ± 0.54. 
From the relation m s {ss) p = — 2(a 2 + a 3 )m s and using m s = 150 MeV, we obtain a 3 m s = 
— (220 ± 40) MeV, which is in the middle of our range of values. 

In order to determine the kaon partition function we generalize the finite temperature 
procedure outlined by Kapusta p4| and Benson, Bernstein and Dodelson [p5[| . First, by 



studying the invariance of the Lagrangian under the transformation — > if ± e ±ia , the 
conserved current density can be identified. The zero component, i.e. the charge density, is 



j = i-fK + d K~ + 2 1 ibK + K- . 
Next, we transform to real fields 0i and 02, 

K± = (0!±i0 2 )/V2. 

The Lagrangian becomes 

C K = \Ad v fad v fa + \Bd v <j) 2 d v <j) 2 + CdvfaVfc 



(7) 



(8) 



where 



.4 



+71 



<t>\ + i4>l 



60i do 02 - s(m# + c)(0i + 
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01 + 70i 
0? + 0| 



C 



1 -7)0102 
0? + 0| 



(9) 



(10) 



Since, to our order of approximation, the counterterm Lagrangian, £ ct , is to be evaluated 
with constant mean fields it plays no role in determining the conjugate momenta. These are 

7Ti = Ado4>\ + Cdo4>2 - 71 602 

7r 2 = Bdo4>2 + Cdo<p\ + 7i&0i • (11) 

2 

The Hamiltonian density is Hk = 7i"i<9o0i + 7T2<9o02 — Ck, and the partition function of the 
grand canonical ensemble can then be written as the functional integral 

Z K = 



[d7n][d7r 2 ] / [d0iP0 2 ] 

J periodic 



X 



exp Udrf d 3 x (^ 1Tl ~^~~ + 



dr 



T-C K (4>i,TTi) +/ijo(0i,7Tj 



;i2) 



Here /i is the chemical potential associated with the conserved charge density, (3 = T _1 is 
the inverse temperature and the fields obey periodic boundary conditions in the imaginary 
time r = it, namely 0j(x, 0) = 0j(x, i3). 



The Gaussian integral over momenta in Eq. fll2[) can be performed yielding 
Z K = N I \dd)i}\dd)ohte- s with 



N Jld^ldfc^e 



dr 



dr dr 



+A(V4> 1 ) 2 + 5(V0 2 ) 2 + 2CV0x ■ V0 2 + 2z( 7 /i + 7 i6)0 2 d T fa 



+ 



71 (m 2 K + c - 2 fib) - 7/i 2 



2C 



at 



(13) 



The dependence of the action ( ]T3|) upon the chemical potential // can also be obtained by 
gauging the fields as discussed in Appendices A and B. 

Next the fields are decomposed into constant condensate parts and fluctuations according 

to 



>i — 



+ # = fd cos a + M £ e^^hAp) 



V 

n,p 

2 = 02 + 2 = fOsin a + e t(p - x+ ^ T) 02,n( P ) • (14) 

V v n,p 

Here V denotes the volume of the system. In the Fourier expansion of the fluctuating part 
0i,o (p = 0) = 02,o (p = 0) = and the Matsubara frequency u n = 2nnT. Since we wish 
to go to one loop order, C ct is evaluated with constant condensate fields and we write the 
contribution as a potential V ct . The remainder of S must be expanded to second order in 
the fluctuations. The partition function can then be written in the form 

Z K = Ne -mv t rM) J n[rf0i,„(p)]M0 2 ,„(p)]7^- 5 — . (15) 
Introducing the definitions 



71, p 



e(0) = 2(m 2 K + c - 2/jb) f 2 sin 2 ±0 - ±/i 2 / 2 sin 2 6 
de 
djj, 



p{9) = -— = Abf 2 sin 2 ±0 + fif 2 sin 2 6 , (16) 



the tree potential is 

Vtree = e(0) • (IT) 

The one loop contribution to the action then becomes 



Sloop = ,-n(-p),02,-n(-p))-D 

n,p \ 1 



>l,nlPj 
>2,n(p) 



where 



(3- 2 f 2 D u = Af{u 2 n + p 2 ) - uj sin 2a + e"cos 2 a + e'^ 

/r 2 / 2 £> 22 = Bf 2 (u 2 n + p 2 ) + u n 9 sin 2« + e" sin 2 a + e'^ 

(3~ 2 f 2 D 12 = Cf 2 (oo 2 n + p 2 ) + 2oo n (J; + 6cos 2 a (j^)' ) + (V - ^ ) sum c-osn 
r 2 / 2 ^2i = Cf 2 (u 2 n + P 2 ) - 2uo n (JL + 9sm 2 a (£)') + {f ~ \* ) sin a cos a , (18) 



with a prime denoting a partial derivative with respect to 9. Here the quantities A, B and 
C of Eq. ( PD are to be evaluated with the condensate fields, giving 

A = 1 + (7 — 1) sin 2 a ; B — 1 + (7 — 1) cos 2 a ; C = (1 — 7) sin a cos a , (19) 

and 7 = (sin^/f?) 2 . The latter replacement should also be made in the measure of the 
integral flT5|), to the order considered. Then the integration can be carried out and the 
grand potential Qk is given by 

2^ = ~ In Z K = V tree + V* + ^ £ In ^det£> . (20) 

Evaluating the determinant yields 

-detD = /3 4 [(w 2 + p 2 ) 2 + 2w 2 (ci + c 2 + c 3 ) + 2p 2 (c 2 + c 3 ) + 4c 2 c 3 
7 L J 

= (3 2 {ul + E 2 + ).(3\u; 2 n + E 2 _), (21) 



where 



£±(p) = P + ci + c 2 + c 3 ± V / c 1 (2p 2 + ci + 2c 2 + 2c 3 ) + (c 2 - c 3 ) 2 . (23) 

Notice that there is no longer any dependence on the condensate phase a. The sum over 
the Matsubara frequencies in Eq. (F2DJ) can then be evaluated in the usual way [2~1 to yield 



VLkV 1 — Vtree + Vet + Vz p + Vthermal with 

KP = /(2^^ [jE+(P)+E - (P)] ^ 

= P~ 1 /-^ HI ~ e-^ +(p) )(l - e-^-W) • (25) 



We are interested in zero temperature here, so Vthermal m &y be dropped. 

The zero-point energy V zp , or equivalently, the one-kaon-loop contribution to the effective 
potential |26| , can be represented as an infinite series of diagrams, some of which are shown 
in Fig. 1. Vertices are given by the Lagrangian (^). External legs represent nucleon or kaon 
mean fields. 

It is interesting to note the chiral order of the diagrams contributing to the effective 
potential. Consider the scattering amplitudes corresponding to the diagrams of Fig. 1, 
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where the external lines are on their respective vacuum mass-shells. The chiral order, v, of 
a given scattering amplitude is given by the Weinberg P7f counting rule, according to which 



an amplitude involving external nucleon lines and external kaon lines is proportional 
to Q u , where Q is a characteristic small momentum scale involved in the process and 

v = 2 + 2L-\E N + Y,{ d i + \ n i- 2 ) » ( 26 ) 

i 

where L is the number of loops, the sum is over all vertices i, d{ the number of derivatives 
that act on the iih vertex and rij the number of nucleon lines attached to the ith vertex. In 
this case Fig. la gives v — 4. Furthermore, it easy to convince oneself that adding a vertex 
from corresponding to the b term in Eq. (||), will decrease v by 1. On the other hand, 
adding a vertex from & 2 \ corresponding to the kinetic, mass, or c terms in Eq. (^) leaves 
v invariant. The fact that these scattering amplitudes give — oo < v < 4 due to the C^- 1 ' 
interaction might appear to be grounds for concern. However, exactly the same apparent 
counting problem arises when one considers the computation of the tree-level in-medium 
kaon propagator. The correct in-medium propagator is indeed obtained by summing an 
infinite series p8j in diagrams that here correspond to amplitudes with decreasing powers 



of v. The evaluation of the effective potential is entirely analogous. We conclude that while 
Eq. (^Bp provides an effective way to organize terms in order of importance for a given 
scattering process, it cannot be naively applied in the same way to the effective potential in 
the case at hand. 



B. Zero-point Divergences 

As shown in Appendices A and B, V c t cancels the ultraviolet divergences in V zp . For our 
purposes the zero-point energy is most conveniently written 

Vzp = ] jllWf (p 2 + Cl + c 2 + c 3 + VV + 2c 2 )(p 2 + 2c 3 )) 2 • (27) 

We regularize the infinities that occur here by introducing an ultraviolet cutoff A and re- 
quiring |p| < A. (Dimensional regularization is discussed in Appendix B.) The counterterm 
potential V ct following from C ct (see Appendix B) is 
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Vet = ~{di + h) - §( Cl + 2c 2 + 2c 3 )(d 2 + / 2 ) + §[(c 2 - c 3 ) 2 + |( Cl + 2c 2 + 2c 3 ) 2 ](d + f ) . 

(28) 

Here we have introduced finite terms f , f 2 and / 4 as well as the quartically, quadratically 
and logarithmically divergent integrals 

f d 3 k , A 4 
rf 4 = / 77, — = Tru (29) 



fc <A (2tt) 3 8tt 
/• d 3 A; 1 A 2 
d2 ~ J k <Aj2^2k- 8^ (30) 
f d 3 k 1 1 , /A\ 



/fe<A (2tt) 3 4A; 3 8tt 2 Vk. 

Here k is an arbitrary scale introduced so that the argument of the logarithm is dimension- 
less. Note that the quadratically and logarithmically divergent pieces of V zp in (|2~?D and 
Vd in (p^ ) are individually dependent upon the chemical potential /i. This dependency on 
the chemical potential \i arises from the functions of the kaon fields multiplying the kinetic 
terms in Eq. (^). However, /i-dependent divergences in the terms V zp and V ct cancel, as 
do the //-independent divergences of the vacuum theory. The cancellation of /^-dependent 
divergences follows from a careful treatment of the counterterm lagrangian which regulates 
the vacuum theory. Further discussion of this point is given in Appendices A and B. Here 
we note that \x may be considered a constant external vector potential A u , with vanishing 
spatial components. The general statement being made is that all Green functions with 
constant external A v fields are finite. It is easy to check this for simple examples related 
to our theory, such as scalar electrodynamics. The fact that it also holds for a complicated 
nonlinear theory expanded about an arbitrary point 0i, 2 in function space is almost cer- 
tainly a consequence of the Ward-Takahashi identities associated with the U(l) symmetry 
transformation preceding Eq. (7). 

In order to determine the finite contributions fi we apply three conditions in the vacuum. 
First Qk, or equivalently the pressure, should be zero in the vacuum and, since there is no 
condensate, the loop contribution must vanish. Secondly, we require the zero-momentum 
propagator to be of the usual form (a; 2 — m^) _1 , see Appendix C. Thirdly we require that the 
four-point vertex (K + K~) 2 /(Qf 2 ) of Eq. (C3) be unmodified in vacuum. These conditions 
imply 
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v zp + v ct 



oe 2 ' 
d 4 



{V zp + V ct ) 



Expanding to fourth order in 9, one finds 





. 



(32) 



u 

fo 



rn 



K 



64vr 2 
1 



rn 



K 



16tt 2 



57T 



In 



m K 



+ c 



(33) 



where using our three momentum cutoff, the constant £ = In 2 — 1/4 = 0.443; interestingly 
a closely similar value for the constant is obtained with dimensional regularization, namely 
( = 1^(3) = 0.461 (Here \I/ represents the digamma function [2~9|). Notice that the arbitary 
scale k only occurs in the combination (d Q + / ) where it cancels out. The disappearance 
of the arbitary scale is one justification for introducing the finite counterterms, we discuss 
another below. 

In evaluating the grand potential for the general case numerical integration of Eq. (|27|) 
for V zp is required. However below threshold when 9 = (or just at threshold when 9 is 
vanishingly small) the expressions simplify greatly and an analytic form can be obtained. 
Explicitly we find C\ = 2(6 + /i) 2 , c 2 = c 3 = \{rn 2 K + c — 2fib — /i 2 ), and eigenvalues 

E± = y/p 2 + «£ 2 ±(b + fi) with £ 2 = b 2 + m\ + c for 9 = . (34) 

Notice that (E± =F /x) is simply the tree level energy of a K ± meson. In this limit the one 
loop energy becomes 

d 3 p 



v zp + v 



ct 



(2tt) 3 
1 



p 2 + e - (d, + h) - i\d 2 + f 2 ) + \e{d, + fo 



64vr 2 



:<-sa«-e 2 )+2e 4 in 



rn 



K 



for 9 = 



(35) 



which vanishes in the vacuum since £ = ttlk there. 



C. Nucleon Hamiltonian 



To describe the nuclear kinetic energy and interactions we employ the non-relativistic 
approach suggested by Prakash, Ainsworth and Lattimer 1 30f1 , which has been rather widely 
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employed, e.g. Refs. [§,0. Here the Hamiltonian density is written 



U N = ~E F u3n + V(u) + n(l - 2x) 2 S(u) , (36) 

where nuclear saturation density is denoted by tiq (no = 0.16fm~ 3 ) and the corresponding 
Fermi momentum and energy are p° F = (37r 2 n /2)3 and E° F = (p° F ) 2 /2M, respectively. The 
potential contribution to the energy density of symmetric nuclear matter is denoted by V(u), 
where u = n/n and it is given by 

V(u) = \Au 2 n, + - - B , ua \ + 3un ^ Q J [-f- - arctan fj . (37) 

Here pf is the Fermi momentum which is related to p F by pf = p° F u^. The remaining 
quantities are parameters which are given in Refs. [|T^,[!tJ. We choose the set which yields 
a compression modulus K = 240 MeV for equilibrium nuclear matter since this is in the 



middle of the generally accepted range of 200 — 300 MeV. The remaining term in Eq. ( |36| ) 
is the symmetry energy, i.e. it describes the change in the energy density when the proton 
fraction x = n p /n differs from |. Again following Ref. |30| we take 



S(u) = (2§ - l) §£° («§ - F{uj) + S F(u) , (38) 

where So — 30 MeV is the bulk symmetry energy parameter. Since neutron star properties 
are not strongly dependent on the form of F(u), we choose the simplest of the suggested 
parametrizations namely F(u) = u. 

In terms of the Hamiltonian density, the grand potential for nucleons is 

^ = Hn - [i n n + xfin , (39) 

where /i = ji n — fi p and the chemical potential /i n is the Fermi energy for neutrons. 



D. Equilibrium Conditions 

We shall study cold, catalyzed neutron stars which are in chemical equilibrium under 
beta decay processes. The process p + e~ «-> n + v e in equilibrium establishes the relation 

[L = fa - fi p = lJL e , (40) 
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since the neutrinos have left the star and their chemical potential is therefore zero. At 
densities where \i exceeds the muon mass m M , muons can be formed by e~ <-> /i~ + + v e , 
hence the muon chemical potential is /x M = \x e = fi. Negatively charged kaons can be formed 
in the process n + e~ <-> n + If - + z/ e when fj, K - becomes equal to the energy of the lowest 
eigenstate of a K~ in matter. Chemical equilibrium requires that /Ik- = A 4 , as we have 
implicitly assumed in the text. 

To complete our grand potential we need the lepton component for which it is sufficient 
to use the standard non-interacting form: 

ih +^H--,){^f(2t 2 + l)tv / ^-ln(t + v / ^)l-4^l (41) 



V 12tt 2 8tt 2 



3tt 2 



where r)(x) is the Heaviside function (r)(x) — 1 if x > 0, 77(3;) = if x < 0) and the quantity 



t = jj? — m^/rrifj,. The total grand potential is then 

fttotai = n K + n N + n L . (42) 

The values of the chemical potential //, the proton fraction x and the condensate amplitude 
6 are determined by extremizing fltotai- 

Extremizing with respect to /i is equivalent to the requirement of charge neutrality. It 
yields 

= -n K +n p -n L , (43) 

where n K and n L are the number densities for negatively charged kaons and leptons respec- 
tively. The kaon charge density is 

n K = f(4bsm 2 \6 + /isin 2 fl) - d{Yzp + Vct) . (44) 

Notice that the fluctuations carry no net charge in the absence of a condensate, i. e. tlk = 
for 9 = 0. Here, and in the other cases, the zero point contribution requires a numerical 
integration and it is straightforward, though tedious, to evaluate the integrand. The lepton 
charge density is 

"* = £ + "<M-^>Rl^- (45) 

Next extremizing with respect to x gives 

15 



= n(2 ai m s - n) sin 2 ±> + ^ +n/I - 4n(l - 2x)S(u) . (46) 

(J Ju 



Finally extremizing with respect to 9 yields 

= f sin 9{m\ + c-2fib-fi 2 cos 9) + Vct) . (47) 

Ou 

Below the threshold for kaon condensation 9 = and this solution is always possible since 
the expression is an odd function of 9. Above threshold the solution 9 > minimizes the 
grand potential. At threshold, where 9 is vanishingly small, it is possible to obtain an 
analytic expression by retaining terms through order 9 2 in Qk- The extremization condition 
becomes 

= (m 2 r + c)(l + 5)-26 y u(l-25)-/i 2 (l + 5) where (48a) 

' "2 2 , t 2i m K 



*=48^(5 -< +fln #J ' (48b) 
This may alternatively be obtained by computing the inverse propagator, D~ 1 (fi), to 
one loop order, see Appendix C. Notice that, by construction, the inverse propagator is 
(fj 2 — m 2 K ) in vacuum since b = c = 5 = there. In fact for any density 5 < 0. That this 
is necessary may be seen as follows. The arguments of the inner square root in Eq. ( p7|) 
is positive for all p, only if C2 = C3 or both C2 > and C3 > if C2 7^ C3. However in our 
calculations we have found that going even slightly above threshold implies C2 7^ C3, so that 
it is the latter condition that must be satisfied. Now the above equation ( |48a| ) can be put 
in the form (1 + 5)(c2 + C3) = — 6bfi5. Since \S\ <C 1 ( Sec. Ill ), and b and fi are positive, 
we must have 5 < 0. It is interesting to note that in the absence of finite counterterms the 
requirement 5 < can only be achieved in a very limited region of parameter space. 

Having satisfied the conditions (f£J), ( [4"6] ) and (pETp, the energy density can be obtained 
from 

r. ^total / Ar\\ 

£ = + /i n n . (49) 

The chemical potential /i n may be obtained by functional differentiation of the energy density 
with respect to the neutron density, viz. \i n = d£/dn n , whence the pressure can be calculated 
from P = ~Vt tota i/V. Equivalently one can write P = n(d£/dn) — £. 
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III. EVALUATION OF KAON LOOP CONTRIBUTION 



In this section, we give numerical estimates for the effect of the loop corrections upon 
the critical density for kaon condensation and upon the properties of neutron star matter 
both above and below the condensate threshold. 

As follows from Eqs. (ESp and (CI), the critical density for kaon condensation is the 
density at which the lepton chemical potential « (an increasing function of density), and the 
position of the zero-momentum K~ pole of the kaon propagator in normal nuclear matter 
(a decreasing function of density) become equal. The chemical potential u and the proton 
fraction x, prior to condensation, are obtained from Eqs. (23) and ([K]) by setting 6 = 0. 
The pole position of the kaon propagator is a solution to the quadratic equation D^iuj) = 0, 
where D" 1 ^) is given in Eq. (C4). 

In Table I, we list the critical density for kaon condensation in units of nuclear matter 
density, /no, for three representative values of the parameter a^m s . In column 2, 

u l cr is the critical density obtained at tree level [13|, corresponding to setting 5 = in Eq. 
(|48|) and omitting the V zp + V ct contributions from Eqs. ([43]) and fl46|) . In column 3, the loop 
contribution is included and the corresponding critical density is denoted by u 1 ^ 1 . We see 
that u 1 ^ 1 < u l cr and, while the magnitude of the threshold shift increases for larger values of 
|a3m s |, it is always small. This due to the fact that in Eq. ( |4"8| ) \5\ < m| c /(487r 2 / 2 ) = 0.06 
which is small in comparison to unity. It is rather easy to show that u^ 1 < w* r . The tree 
level inverse propagator (/}*(«) ) _1 = — (c2 + C3) = at threshold, while with loops included 
one has —(1 + S)(c2 + C3) = 66u5. Since D t (fi) is a monotonic function of density and we 
have pointed out that S < 0, it follows that u 1 ^ 1 < w* r . 

The change in the threshold due to loops, small though it is, is dominated by the 
term. This is shown in column 4 of Table I for which we set c = in the loop expressions. 
In this case the threshold is essentially unchanged from tree level. The effect of the c term 
can become quite critical at densities above threshold and in Eq. (||) the scalar density 
n s = (NN) has been replaced by the baryon density n = (N^N). In Ref. [ffl it was argued 
that in the heavy baryon limit, and to the chiral order considered, no distinction should be 
made between the operators NN and N^N in the Lagrangian ([I]). Corrections to this would, 
in principle, require other effects of similar order for consistency. While n s ~ n for n ~ n Q , 
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at high density n s becomes less than n. The results of Ref. |3T| indicate that this can be 
significant, in particular it leads to an increase in u cr . Since our results are sensitive to this 
question, in some of our calculations we shall replace the baryon density n in the quantity 
c of Eq. (H) by a phenomeno logical form of n s . Based on a simple parameterization of the 



effective nucleon mass obtained by Serot and Walecka [3~2~|, we use 

2.5213 / 1 



n, = n- 



u V l + exp( Ml/3 - Lm ^ 



(50) 

0.4417 ) J 

for u > 1.0342 and set n s = n for u < 1.0342. This is extremely crude, but it will give some 
qualitative feel for the effect. The results for the critical density with this approach, labelled 
c ns , are shown in columns 5 and 6 of Table I. The shift in the critical density between tree 
and tree + loop levels is slightly smaller, but basically similar to that obtained before. The 
overall increase in the critical density is expected, but is somewhat larger than models which 
treat nucleons relativistically would suggest. 

Even below threshold there is a kaon zero-point contribution for normal neutron star 
matter. Firstly there is a shift in the minimum due to the zero point contribution to the 
equilibrium Eq. (||j) at 9 = 0. Secondly there is the zero point energy (^) itself. The 
shift in the proton fraction x and the lepton chemical potential \i is negligible. Table II 
shows that e zp , the zero-point energy density is approximately 1-3% of the total. Thus the 
matter properties are only affected to a small degree so that it is not necessary to refit the 
parameters of the nuclear Hamiltonian. 



Above the critical density the properties of matter are obtained by solving Eqs. (f!3"D, 
(JlH), and fl47f ) numerically. In practice, one uses a three-momentum cutoff A to evaluate 
Eq. fl2"7|) and derivatives thereof. The cutoff used was A = 3 x 10 4 MeV/c, and results were 
found to be extremely stable with respect to variation in A. For example taking A = 1 x 10 4 
MeV/c, the largest change in any of the tabulated values was 0.3%. 

In Table III, we list the properties of normal neutron star matter (column 2) and matter 
containing a kaon condensate at tree level (columns 3-7) and at one-loop level (columns 
8-13) for a^m s = —134 MeV. Here the kaon fraction is defined as xk = nx/n in terms of 



the kaon density rix of Eq. (44), and the lepton fraction xl = ni/n may be obtained from 
xl = x — Xk- The energy density is also listed. We have not felt it worthwhile to tabulate 
the remaining thermodynamic variable, the pressure, in view of the small size of the loop 



effects. For u > 9.03 the parameter C3 becomes negative which precludes real solutions to the 
zero point energies; we discuss this further below. In Table IV, we give the corresponding 
quantities for 03/77.5 = —222 MeV, in which case C3 becomes negative for u > 5.07. This 
difficulty does not arise if we employ the phenomenological scaling given by Eq. (|50| ) and 
the results are given for d^rris = —222 MeV in Table V. If one simply sets c = in the 
loop calculations the results shown in Fig. 2 are obtained. In this figure Ae represents the 
energy density gained by the formation of a condensate and is the difference between the 
total energy density and the tree level energy density for normal nuclear matter. 

We first note that loop effects for the case a^m s = —134 MeV are essentially negligible. 
Although they increase in magnitude for larger values of |a3m s |, only small modifications of 
the tree level results are found. Table IV shows that just above threshold the tree + loop 
calculation gives values of 9 and x which are greater than, and fi which are less than, the 
corresponding tree level results. This naturally follows from the fact that u^ 1 is less than 
u' r since, once a condensate is present, 6 rises rapidly and x is enhanced, while /1 is reduced 
At high densities, this ordering is reversed. Near the limiting density, u ~ 5.07, 9 typically 
decreases by ~ 3% while x decreases somewhat less, ~ 1%, and \i increases ~ 6%. The 
zero-point energy density is of the order —(1.2 — 3.5) MeVfm -3 and can decrease the tree 
level result significantly at high density when the latter can be small. For a given density 
the effects are smaller if the scalar density of Eq. (^) is employed as shown in Table V. 
On the other hand, setting c = altogether as in Fig. 2 yields larger effects at very high 
density. For a 3 m s = —310 MeV, the above trends are also present, and are enhanced, as 
expected. The coefficient C3 goes to zero for u > 3.32, at which point e zp ~ — 7.2MeVfm~ 3 
and e ~ -8.5 MeV fm" 3 . 

As mentioned above, C3 often vanishes at some density above the critical density. The 
condition C3 = implies de/d9 = 0, where e is the tree level energy density function of 
Eq. (|T7]). The significance of this condition is best illustrated by the simple 0(2) invariant 
Lagrangian C = \{d v (pi) 2 + \{d v $2) 2 — U(<fi), where (ft 2 = (p\ + $\. The one-loop effective 
potential [|26| is U l (<f)) = zTrln((A; 2 + U")(k 2 + U'/(f))), where the trace is over momentum. 
Now take U(4>) = \m 2 4> 2 + A0 4 , with A > 0. For m 2 > 0, U l is real for all 0. On the other 
hand, for m 2 < 0, the "Mexican hat" potential, the one-loop effective potential is complex in 
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the region U' < which signals a physical instability. Weinberg and Wu [|33| show that the 
imaginary part can be interpreted as a decay rate per unit volume of a well-defined state. 
The connection between the instability and the chaotic dynamics of a classical field theory 
has recently been discussed by Matinyan and Miiller ||34|| , see also references therein. 

Unfortunately no solution is known to this general problem. Our case is entirely analo- 
gous. At threshold = C3 > corresponds to m 2 > 0. As the density is increased, the tree 
level potential eventually corresponds to m? < 0. This effect is more pronounced for larger 
values of |a3m s |. Since the effective potential is by definition real, this shows that some, as 
yet unknown, technique is needed in these regions of parameter space. The two other cases 
we have considered, excluding c from the loop contribution and employing the softer density 
dependence of c with the parametrized scalar density n s alleviate this problem. 



IV. CONCLUSION 



We have derived the zero-point energy due to fluctuations in kaonic modes in neutron 
star matter containing a kaon condensate. To our knowledge, this calculation is the first 
in which the nonlinearities of the chiral Lagrangian have been taken into account at the 
one loop level for densities above the critical density for kaon condensation. The zero-point 
energy was derived through the thermodynamic partition function by considering fluctua- 
tions around an arbitrary condensate. We have demonstrated by gauging the counterterms 
of the vacuum theory that no divergences are present when matter is introduced via the 
charge chemical potential. In the usual way the total thermodynamic potential, including 
the classical contribution, is extremized so that the minimum energy ground state is found 
for locally charge neutral and beta-equilibrated matter. 

The critical density for kaon condensation was found to be reduced by fluctuations, 
but only by less than 1% for the range of parameters considered. At higher densities the 
loop effects were negligible for a 3 m s = —134 MeV which corresponds to zero strangeness 
content for the nucleon. The effects were slightly increased with larger values of |a 3 m s |, 
corresponding to larger value of the kaon nucleon sigma term E^^, but they remained small 
in comparison to the dominant tree level contributions. The size of the loop effects relative 
to the tree contributions is set, at least in the threshold region, by the parameter 5 which is 
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less than m 2 K / (48ir 2 f 2 ) = 0.06. In some cases it was not possible to obtain solutions up to 
arbitrary density because the effective potential became complex- this is a well-known and 
unsolved problem in one loop calculations with spontaneous symmetry breaking. 

The fact that changes induced by the zero-point motion are generally small leads one 
first to conclude that no major revision is needed of earlier investigations in which such 
fluctuations have not been taken into account. Second, the loop expansion appears to be 
convergent where it is well defined. The formalism presented here yielded the contribution to 
the grand potential arising from thermal excitations of the fluctuation modes. These effects 
would have to be included at finite temperatures and would be involved in the extremization 
conditions. Such an analysis is needed for a consistent study of newborn neutron stars. 
Extension of the present work to include the full octet of pseudoscalar fluctuations, and to 
include additional terms in the chiral expansion could also be contemplated. 
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APPENDIX A: CANCELLATION OF CHEMICAL POTENTIAL 
DIVERGENCES IN A TOY MODEL 

In this Appendix, we illustrate how divergences dependent on the charge chemical po- 
tential \i are cancelled in a toy model which contains two real fields, 0iand 02- The model 
we consider is 



and 7 is an arbitrary analytic function of 4> 2 = {4>\ + 4>\) with 7(0) = 1. It is useful to discuss 
this model to one loop order since it has a simpler structure than the chiral Lagrangian, but 
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ct 



where 




(Al) 
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the /i-dependent divergences are removed in a similar way. We first derive C c t, and then 
demonstrate that on introducing an external chemical potential //, the effective potential 
(or equivalently, the zero-point energy) is free of divergences. Finite counterterms are not 
relevant in the present context and will be ignored. 

Expanding the function 7 to all orders in the fields gives the the standard kinetic energy 
term plus an infinite set of interaction terms involving an even number of fields and two 
derivatives. From these interaction vertices, we may compute one-particle irreducible 2n- 
point Green functions to one-loop order. The counterterm Lagrangian C ct is defined in such 
a way as to minimally cancel divergences. 

Consider first quadratic divergences, proportional to d 2 (defined in Eq. (p0[)). The 4- 
point vertex contributes a quadratic divergence to the 2-point function proportional to d 2 
and the external momenta. The corresponding counterterm Lagrangian is proportional to 
d2[{dy(j)i) 2 + {d v (f)2) 2 }- The 6-point vertex contributes to the 4-point function so as to require a 
counterterm Lagrangian proportional to d 2 <f) 2 [(d l j<f)i) 2 + (<9j,02) 2 ]. Continuing this sequence, 
we find that the contribution of quadratic divergences to the counterterm Lagrangian is 
—^2/(47) • (l" + 1' / 4>)[(d u (pi) 2 + (<9i/02) 2 ], where primes denote derivatives with respect to 0. 

Now consider logarithmic divergences, proportional to do (see Eq. Q3"TP). Green functions 
containing such divergences typically contain two vertices, and evaluating them requires 
involved combinatorics. It is simpler instead to derive the counterterm Lagrangian using 
the effective potential (see e.g. Ref. PEfl). One thus avoids combinatorics, and the method 
is more easily applied to Lagrangians more complicated than Eq. (|A1|) . 



where 0i and 02 are arbitrary fields which, for the present, may vary in spacetime. Expand- 
ing to quadratic order in the fields 0' x and 0' 2 leads to 



Let 



01 + 01 , 02 -> 02 + 2 



'2 > 



(A2) 




(A3) 



where the matrix of second order variations is given by 



D = -jOl+L-d-M , 



(A4) 
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L v = 2 



M 



7 



02*^01 02^02 



(^0O 2 + {d v fa 



<Pll" + <P% 0i02(7" -J) 
0102(7" -i) 0h" + 0?i 



(A5) 



(A6) 



Integrating out the fields 0' x and 0' 2 gives the unrenormalized effective potential to one-loop 
order 



— -Tr \nD = — zTr In 7 — zTrln(— □) 



2&k \-2ryn 



[0 -L-L-d +2M+(d-L)} 



(A7) 



where, as before, the trace is taken both over spacetime and charge. The factor — iTrln7 is 
cancelled by a corresponding contribution from the measure, and will not be considered fur- 
ther. In the last expression, we have symmetrized the derivative contribution and separated 
the mass-like contributions involving [2M + (d ■ L)]. We shall see that only the mass-like 
contributions are important for our discussion. 

The trace over spacetime is evaluated by inserting complete sets of four-momentum 
eigenstates. A momentum expansion is performed on the resulting integrals, leading to 
a factorization into ultraviolet divergent loop integrals and integrals over functions of the 
fields. To minimally cancel the divergences in the effective potential we need the following 
counterterm Lagrangian 



p pquartic . pmass—like , ^derivative 



(A8) 



where the contribution from the quartic divergence (|29|) is 

J d 4 x £T rUc = -iTrln(-D) = J d 4 x d 4 



(A9) 



The contribution from mass-like terms is 

2 



<TxC 



mass— like 
ct 



2 fc=l k 
47 



[2M + (d ■ L)] 



[ d 4 x Tr[2M + (d ■ L)\ - / d 4 x Tr[2M + (d ■ L)\ 

J ' luT J 



I67 2 



-^Jd'xU 1 '' + 1 -\ [{d^f + (DM 2 
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dn 



4 7 2 

i 



d x 



(dM 2 + (d v <f> 2 ) 



7 



+ 



7 



(dM 2 + {d^f 



7"(01^0i + 02^02 



7' 



(dyfadxfa + ^0 2 9 A 2 ) 2 + (01^01 + 



+2((<9 i ,0 1 ) 2 1 D0 1 + {VfadMfaUfa + {1^2} 



(A10) 



For the contribution from derivatives one finds (see also Refs. [35,36]) 

j d , x privative = * ^ ^ (j-(L- d~d -L)) 

= ~J d*x Tr^(d u L x - d\L u ) 2 - ^(d u L x - d x L v )[L x , L v \ 

+^L 2 V L 2 + ^(^a) 2 } ~^Jd A x TrLl . (All) 

This concludes the derivation of C c t, to one-loop order. 

Having derived a one-loop finite vacuum theory, we are now in a position to study the 
theory in the presence of a finite chemical potential /i, at finite inverse temperature j3 = 1/T. 
As before, the Euclidean time is r = it, and the action e tS —>■ e~ s under the transformation 
to Euclidean space. The chemical potential acts as an external vector field with spatial 
components zero. It enters by "gauging" the fields 0i and 02 via the transformation 



■ d(pi .901 
i— > i— /U02 ; 



dr dr 
The resulting action is given by 



• <902 .d(f>2 , , 
l-^ > l-^ h /Mpl ■ 

or or 



(A12) 



S = S-y + Set , where 



S, 



. — | / dr I d 3 x ^C^ . 



—dd — 



d 2 
4 7 



Set — I dr I d 3 x C ct and 



(7'T + 



I67 2 



7 



(AJ 2 + ■ ■ 



c„ = (9 r 0x) 2 + (9 r 2 ) 2 + (V0O 2 + (V0 2 ) 2 + 2^(0 2 a r 0! - 0x^02) - m 



2 J.2 



(A13) 
(A14) 

(A15) 
(A16) 
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Eq. ([A lip is invariant under the transformation (|A12 ) and is therefore not included above. 



We have now derived the action to one-loop order, for arbitrary fields 0i and 02 in the 
presence of a chemical potential \i. 

We then take the final step, and evaluate the one-loop effective potential for fields 0i 
and 02 which are constant in spacetime (see Eq. (0)) so that their derivatives vanish. 



Expanding Eq. (|A13|) to quadratic order in the fields and integrating out the fluctuating 



field gives the action 



S — Stree ~\~ Sl oop + S c t , 

P 

Stree = ~\ I J d 3 X -ffi 2 ^ , (A17) 


S loop = -llndetD , (A18) 
where the matrix of second order variations about the constant fields, D, is given by 
\D U = 7 (^ 2 + p 2) _ 2u n yL(j> l( j ) X - /i 2 7 - \Act>\ + 40?) - ±/i 2 2 7 " 



P 

j 2 D 22 = 7 (^ 2 + p 2 ) + 2u; n //0 1 2 ^ - /i 2 7 - \v 2l -(<Pl + 402) - f/^W 



—D 21 = -2fiu n 7 + (fiJ- - V0X02 / + 3^- , (A19) 



j 2 D l2 = 2fiu n (7 + 2 ^ - l / i 2 1 2 (7" + 3 f\ 

IP' y 

and u; n is the Matsubara frequency. The determinant is given by Eq. (|2l]) , with an additional 
factor of 7 _1 on the left, and with 

ci = —(7 + 07') , c 2 = -£-(27 + 4 7 '0 + 0V) , c 3 = -^(2 7 + 7 '0) . (A20) 
7 47 47 

As before, it is convenient to factorize the determinant, leading to dispersion relations for 

two modes given by Eq. fl23l). Eq. ( A18 ) then separates into a zero-point and a thermal 



contribution 

P 



Sloop = J dr J d 3 X (V zp + Vthermal) , (A21) 



with V zp and Vthermal given by Eqs. fl24l ) and ([250 respectively. The counterterm action 
obtained from Eqs. (|A9p, ( |A10| ) and (|A1 1| ) after gauging and setting 0j to be constant can 
be written 
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S ct = fdrf d 3 x Vet , (A22) 
o 

where V c t is given by the divergent parts of Eq. (|28|) , i.e. the terms involving do, d 2 and d^. 
As noted in Sec. IIB, the \i dependent divergences in V zp are exactly cancelled by those in 

Vct. 



APPENDIX B: CANCELLATION OF CHEMICAL POTENTIAL 
DIVERGENCES IN THE CHIRAL LAGRANGIAN 



In this appendix, we give details of the cancellation of quadratic divergences dependent 
on the chemical potential /x for the chiral Lagrangian, Eq. (Q). (The cancellation of quartic 
divergences needs no discussion since it is straightforward.) We follow the procedure given 
in Appendix A. It is first convenient to break the Lagrangian down as follows: 

£-K = C'H + £-v + £>u + £>ct 

C H = \A (d u <Pi? + \ B (d u <j) 2 ) 2 + C dyfadvfa 



C v = -2f 2 (m 2 K + c) sin 2 ±6 = -U(4>) 



(Bl) 



Note that in addition to having more terms than the toy model of Appendix A, the kinetic 
term, Ch, is nondiagonal here. 

Expanding as in Eq. (|A~2j), the matrix D of second order variations of the fields is 
given by 

I A C\ 

D = —TO +L-8-M ; T= . (B2) 

\C B) 

The matrix L can be written 



L 



J v 2 



,7' 



where 



Lh 2 




L v 2 



5 v0 2b^- 




(B3) 
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The matrix M takes the form 



M, 



M = M Hl + M H2 + M Vl + M V2 + M v where 

! r . / 

Hi = -7ri9 (dM 2 + (d„<p 2 ) 2 



M, 



1-7 



Hi 



^□0! 0^02 
^□01 02 D 02 



Mi 



Vi 



7i ' 
2b^ 



-01<9°0 2 -02#°02 



Af % =-6(0l5°02)— 2*2 a * 2 

^ I 0102 (7? - hi) 0M + 



Mr; = — 



2 ^2 

i f <ftu« + ft% hMU"-f, 

^ Ui0^"-f) 0^" + 0?f 



(B4) 



As in Appendix A primes denote derivatives with respect to 0. It may be checked that by 
gauging D according to Eq. (|A12 ) and then taking the fields 0, to be constant, one recovers 
Eq. Q18[). 

It is now straightforward, but tedious, to compute the the quadratically divergent con- 
tributions to the counterterm Lagrangian by evaluating the terms of the form 



iTrtT^M) - ±Tr(T- l L 



1 tv\2 



(B5) 



Omitting terms that are invariant under the transformation (|A12 ), we find the following 
contributions to £ rf /d 2 : 



\Tr{T- l M Hl ) 
iTr^M^) 

\Ti{T~ l M Vl ) 

lTi(T- l M V2 ) 



-\ [{d v <t> x ) 2 + (<^0 2 ) 2 
1-7 



v 7 + 70, 



20 2 



-(0lQ01 + 2 n0 2 ) 



7i 

(0i a 



7i 



i 6 y; + _t ( 0, a u 2 ) 



70 



(B6) 
(B7) 

(B8) 

(B9) 

(BIO) 
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— §Tr(T -1 Ly) 2 = + Hi) (Bll) 

-iTr^- 1 ^ + L^]) 2 = ~i(j- W<h + <hP<h) a (B12) 

-iTr(T-iL» v T-'[L" Hl + =^iX~ ^r) ^° <fc) ■ (B13) 



Following Appendix A, we use the transformation ( A12 ) on C ct and find that the quadrat- 
ically divergent contribution to V ct is 



VT ad = d 2 { \fcf ( 7" + + |/i 2 (l - 7) + ¥<M + I tJ + -J 



+6/3(1 _ 7 _ 0T ') _ + 07i ) — | (V + ^)\. (B14) 

7 722^ 7 0y j 

It may be checked that this agrees with the divergent parts of Eq. (|28| ) using the explicit 
form of the coefficients c\, c 2 and C3: 

Cl = 2/i 2 ( 7 + 7 ') + 2/i6 (271 + + 7 '^ + 26 2 ^( 7 i + M) (B15) 

V 2 7 2/ 72 2 

c 2 = -±// 2 (2 7 + 40 7 ' + 2 7 ") - ^6(271 + 40 7 1 + 2 7 1') + \U" (B16) 



2 

1,. ., V 



c 3 = -}„' (2 + " i/*-(271 + *rj) + ^ ■ (BIT) 

We now discuss the chiral ordering of the above counterterms in the Lagrangian (|TJ). 
To do so, it helpful to evaluate the divergent integrals of Eqs. (p9|), fl3~0|) and (|3l|) using 
dimensional regularization. This gives 

rfl ./|^^._^(I +lD .i_ + | * ( 3)) 
(2^) 3 2 v /A; 2 + m 2 Y Ve m K 2 



' (2ix) 6 A(k 2 + m l K yi l 87r 2 Ve z ' 



Here, we have included explicitly in the kaon propagator. A factor of v47r has been 
subsumed in k which is introduced so that the dimension of the integral is independent of 
the dimensionality of spacetime. The parameter e = 3 — n, where n is the dimensionality of 



the integral, and ^(n) is the digamma function [29|. In minimal subtraction the integrals 
of Eq. ( pi8|) are regulated by removing the 1/e terms. It is then evident that in the chiral 
expansion d±, d 2 and do are 0(q 4 ), 0(q 2 ) and O(q ) respectively. 
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To evaluate the chiral order of counterterms, we promote the nucleoli fields in b to 
dynamical fields. We may then determine the chiral order of counterterms from Weinberg's 
counting rule (p6|), either directly from the tree- level counterterms or from the representation 

_____ ( /Q\ 

of these terms as one- loop diagrams. Counterterms (__f), ( |B9|) and (|B13|) contribute to C ct . 
Expanding to lowest order in the fields <pi and 02, these terms renormalize the kaon-nucleon 

,(4) 
-"ct • 



vector interaction coupling constant. Counterterms (pq ), (__]) and (|B12| ) contribute to C 
To lowest order in the fields, they contribute to wave function renormalization. Eq. ( |B10| ) 
is also 0(q A ), but to lowest order in the fields gives an overall irrelevant constant and a 

,(2) 



mass renormalization. Eq. ( |B11|) contributes to C ct and yields a four-nucleon interaction 
to lowest order in the fields. 

In this Appendix, we have not discussed logarithmic divergences. We believe the cancel- 
lation of n dependent logarithmic divergences will proceed as in the toy model of Appendix 
A. However we have not derived the full counterterm Lagrangian proportional to d since 
this, while in principle straightforward, would involve prohibitively complicated algebra. 
The counterterms will have a chiral order which ranges from to 4 as can be seen from 
Eqs. ( |A10| ) and (|A11|) in conjunction with Eqs. ( P3|) and (__]). In conclusion we note that 
although we have referred to chiral counting throughout this manuscript, our Lagrangian 
should be regarded as a "U(l)" or "0(2)" analogue of the chiral expansion. 



APPENDIX C: CRITICAL DENSITY FROM THE KAON PROPAGATOR AND 
FROM EXPANSION OF ZERO-POINT ENERGY 

In this appendix, we derive the critical density from the kaon 2-point function to one-loop 
order and show that it is equal to that obtained from expanding the zero point energy. 

Close to threshold, the grand potential may be expanded in terms of the condensate as 
follows 

^ = ^-\D-\^)fW + 0^) (CI) 

where D is the zero-momentum kaon propagator, relevant for s-wave condensation, and 
H is the lepton chemical potential prior to kaon condensation. The threshold for kaon 
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condensation is therefore given by D~ l ([i) = 0, i.e. the density at which the energy of 
negatively charged kaons equals the lepton chemical potential. 

From the chiral Lagrangian, Eq. (d), we may read off the interactions. The 2-point 
interactions are 

C U2 = -( m 2 K + c)K + K- , C V2 = -ibK + d° K~ . (C2) 

Here we treat the standard mass term as an interaction and take the kaon propagator to be 
massless. The 4-point interactions are [|[] 

C m = -^- 2 {K + d v K~f , C V4 = ~^- 2 {K + d° K-)K + K~ 



Cu±= { -^fpt {K + K - ) 2 . (C3) 



i 

Qf 2 

We compute the kaon 2-point function with external four- momentum p = (u,0). The 
interactions (|C2 ) contribute at tree-level. To one-loop order, we have the diagrams of Fig. 
3. The contributions of Figs. 3(a), 3(b) and 3(c) are -u 2 /(3f 2 ), -2buj/(3f 2 ) and 2(m 2 K + 
c)/(3f 2 ), respectively. In the last two diagrams, the internal lines contain the interactions 
given in Eq. (|C2j ), which are denoted by m# and b, respectively. Fig. 3(d) contributes 
~(, m K + c )/(3/ 2 ); and Fig. 3(e) contributes 2bu/f 2 . In addition, we must include all 
possible insertions of the interactions ( |C2j ) on internal lines. Including the counterterm 
contribution, the result is 

D~ x {lo) = uj 2 + 2bu — m 2 K — c + (co> 2 — 4ub — m 2 K — c)S where (C4) 

* = -ip (I - 2( * + h) + 2f2(<i ° + A) ) (C5) 

Eq. (|C5|) agrees with the expansion of V zp + V ct given in Eq. (H), as required, and in vacuum 
we obtain the free inverse propagator D~ 1 (oj) = (d 2 — m 2 K ). 

Note that Figs. 3(a), (c), and the contribution to (d) proportional to m 2 K , do not involve 
nucleon fields and are mass and wavefunction renormalization diagrams in the vacuum. 
These diagrams are 0(q A ). If we promote the nucleon fields to dynamical ones, we generate 
nucleon legs on Figs. 3(b) and (e), corresponding to Figs. 1(g) and 1(f) respectively, and 
find that they are 0(q 3 ). These diagrams correspond to Figs. 1(a) and 1(b) of Ref. 0] 
respectively; since only nucleon and kaon fields are used in the present calculation, the 
remaining diagrams of Ref. Q do not enter. 
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FIGURES 




FIG. 1. Some contributions to the zero-point energy (effective potential to one-loop order). 
The dashed (solid) lines denote kaons (nucleons) and the external legs are mean fields. There are 
5 additional diagrams with 6 external legs. 
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u 



FIG. 2. Ground state neutron star properties for matter containing a kaon condensate with 
CL3m s = —222 MeV and c = in the loop contribution. The dashed lines give the tree-level results 
and the solid lines show the effect of also including the loop contribution. Plotted are the charge 
chemical potential fi, the proton fraction x, the condensate amplitude 6, and the energy gained in 
forming a condensate Ae, as a function of u = n/riQ, the density in units of equilibrium nuclear 
matter density. 
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FIG. 3. Contributions to the kaon propagator to one-loop order. The dashed line denotes a 
free kaon propagator. See Appendix C for further notation. 
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TABLES 

TABLE I. Values of the critical density for kaon condensation, u cr = n cr /no, in units of 
equilibrium nuclear matter density, for various choices of the parameter a^m s . Superscript t denotes 
the tree level contribution, and I denotes the loop contribution. The scalar kaon-nucleon interaction 
c is defined in Eq. (5), while c„ s is the corresponding quantity with the density n replaced by the 
scalar density n s given in Eq.(50) 









c = in / 






a 3 m s 




u t+i 


u t+i 






-134 MeV 


4.181 


4.175 


4.179 




5.069 5.069 


-222 MeV 


3.082 


3.064 


3.081 




4.231 4.224 


-310 MeV 


2.417 


2.397 


2.417 




3.473 3.454 


TABLE II. Properties of normal, n,p, e ,/i , 


neutron star 


matter in the presence of kaon 


zero-point fluctuations as a 


function of density expressed in units of equilibrium nuclear matter 


density, u = 


n/riQ. Here a 3 m s = —222 MeV. Tabulated are the proton 


fraction x, the charge 


chemical potential /U (in MeV), the zero-point contribution to the 


energy density e zp (MeVfm 3 ), 


and the total 


energy density 


e (MeVfm" 3 ). 








u 


X 


I 1 




e zp 


e 


0.50 


0.015 


64.8 




-0.01 


1.31 


1.00 


0.039 


110.9 




-0.05 


4.56 


1.50 


0.074 


145.2 




-0.16 


9.28 


2.00 


0.106 


172.7 




-0.37 


15.10 


2.50 


0.132 


195.9 




-0.70 


21.80 


3.00 


0.154 


216.2 




-1.17 


29.23 


3.06 


0.156 


218.6 




-1.24 


30.23 
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TABLE III. Ground state neutron star matter properties for normal n,p,e,fi matter, and 
matter containing a kaon condensate to tree level and one-loop level. Here a^m s = —134 MeV. 
Tabulated are the density in units of equilibrium nuclear matter density, u, the amplitude of the 
condensate 9 (in degrees), the charge chemical potential fi (MeV), the proton fraction x, the kaon 
fraction xk, and the total energy density e (MeV fm~ 3 ). The zero-point contribution to the energy 
density is e zp (MeVfm -3 ). 



n,p,e, fx n,p,e, fi, if (tree) n,p,e,fi, K (tree+loop) 



u 


e 


i_i 
U 


A* 




X 


XK 






II 

U 


A* 




X 


X K 




e 


4.18 


51.98 


0.00 


256.3 


0. 


.194 


0.000 


51 


.98 


0.00 


256.1 


0. 


.194 


0.000 


-0.10 


51.88 


4.20 


52.36 


5.67 


255.3 


0. 


.198 


0.007 


52. 


.36 


6.56 


254.7 


0, 


199 


0.010 


-0.11 


52.25 


4.70 


62.60 


29.55 


226.9 


0. 


.281 


0.166 


59. 


.87 


29.76 


225.4 


0. 


.283 


0.170 


-0.23 


59.64 


5.20 


73.40 


41.02 


197.7 


0. 


.346 


0.280 


63 


.31 


41.01 


195.8 


0. 


.348 


0.284 


-0.35 


62.95 


5.70 


84.74 


49.16 


168.8 


0. 


395 


0.361 


63 


.33 


48.99 


166.8 


0. 


396 


0.363 


-0.41 


62.91 


6.20 


96.57 


55.18 


141.3 


0. 


.431 


0.415 


60. 


.67 


54.90 


139.5 


0. 


.431 


0.416 


-0.40 


60.27 


6.70 


108.87 


59.63 


115.9 


0. 


.458 


0.451 


56. 


.03 


59.31 


114.5 


0. 


.458 


0.451 


-0.32 


55.70 


7.20 


121.61 


62.91 


92.9 


0. 


.477 


0.474 


49. 


.97 


62.61 


92.1 


0. 


.477 


0.474 


-0.21 


49.75 


7.70 


134.78 


65.42 


72.4 


0. 


.492 


0.490 


42. 


.90 


65.20 


72.0 


0. 


.491 


0.490 


-0.11 


42.79 


8.20 


148.35 


67.38 


54.2 


0. 


.502 


0.502 


35. 


.11 


67.25 


54.1 


0. 


.502 


0.502 


-0.04 


35.08 


8.70 


162.30 


68.92 


38.0 


0. 


.511 


0.510 


26. 


.80 


68.87 


38.0 


0. 


510 


0.510 


-0.01 


26.80 
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TABLE IV. Ground state neutron star matter properties for normal n,p,e,fi matter, and 
matter containing a kaon condensate to tree level and one-loop level. Here a?,m s = —222 MeV (See 
Table III for notation). 



n,p,e, \i n,p,e, fi, if (tree) n, p, e, fi, K (tree+loop) 



u 


( 




e 






X 






e 


it 




X 








3.06 


31 


,47 














0.00 


218.6 





.156 


0.000 


-1. 


.24 


30.23 


3.08 


31 


.77 


0.00 


219.2 





.157 


0.000 


31.77 


7.20 


216.7 


0. 


.163 


0.012 


-1. 


,27 


30.50 


3.26 


34 


.82 


22.86 


200.0 





,222 


0.113 


34.10 


24.31 


196.3 





,230 


0.128 


-1. 


.67 


32.44 


3.46 


38 


.35 


33.81 


177.1 


o 


,287 


0.220 


35.13 


34.88 


172.2 


0. 


.296 


0.235 


-2. 


,17 


32.97 


3.66 


41. 


,99 


42.29 


152.8 





.345 


0.308 


34.52 


43.03 


147.0 





.353 


0.321 


-2. 


.70 


31.82 


3.86 


45. 


,75 


49.41 


127.6 





.395 


0.378 


32.36 


49.71 


121.3 





.401 


0.387 


-3 


.18 


29.17 


4.06 


49. 


.60 


55.40 


102.0 


0. 


.438 


0.430 


28.81 


55.15 


96.0 


0. 


.441 


0.435 


-3. 


.56 


25.22 


4.26 


53 


,56 


60.50 


76.8 





.473 


0.470 


24.06 


59.68 


71.7 





.473 


0.470 


-3 


.80 


20.21 


4.46 


57. 


.61 


64.91 


52.6 





.501 


0.500 


18.30 


63.52 


48.9 


0. 


.498 


0.498 


-3 


.90 


14.35 


4.66 


61 


.76 


68.66 


29.9 





,524 


0.524 


11.69 


66.77 


27.8 





.519 


0.519 


-3 


,85 


7.80 


4.86 


66 


,00 


71.80 


9.0 





.542 


0.542 


4.41 


69.53 


8.3 





,535 


0.535 


-3 


.69 


0.71 


5.06 


70. 


.32 


74.43 


-10.3 





.556 


0.556 


-3.40 


71.89 


-9.7 





.548 


0.548 


-3 


.44 


-6.82 
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TABLE V. Ground state neutron star matter properties for normal n, p, e, fx matter, and matter 
containing a kaon condensate to tree level and one-loop level. Here a^m s = —222 MeV and the 
scalar density of Eq. (50) is employed in the scalar kaon-nucleon interaction c„ s . (See Table III 
for notation.) 



n,p,e,fi n,p,e, fi, K (tree) n,p,e,n,K(tiee+loop) 



u 




9 


11 






X 




( 




e 


LI 






X 




z p 




4.22 


52.84 


















0.00 


257. 


.6 





.195 


0.000 


-0.09 


52.75 


4.23 


52.96 


0.00 


257. 


.8 





.196 


0.000 


52 


.96 


2.56 


257. 


4 





.196 


0.001 


-0.09 


52.88 


4.30 


54.36 


8.93 


255. 


.5 





,205 


0.017 


54 


.28 


9.57 


255. 








.206 


0.020 


-0.09 


54.24 


4.80 


64.72 


25.67 


239 


.0 


o 


.263 


0.128 


62 


.86 


25.86 


238, 


.1 


0. 


.264 


0.131 


-0.11 


62.79 


5.30 


75.63 


34.53 


223 


.2 





.308 


0.210 


69 


,54 


34.51 


222, 


.2 


0. 


.309 


0.213 


-0.13 


69.44 


5.80 


87.07 


40.88 


208. 


,2 





.342 


0.272 


74 


.69 


40.73 


207. 


3 





.343 


0.274 


-0.12 


74.59 


6.30 


98.99 


45.75 


194. 


.3 





.369 


0.318 


78 


.63 


45.51 


193. 


5 


0. 


.370 


0.319 


-0.10 


78.55 


6.80 


111.38 


49.58 


181. 


.3 





,390 


0.353 


81 


.64 


49.29 


180. 


.8 





.391 


0.354 


-0.08 


81.58 


7.30 


124.21 


52.64 


169 


.5 





.407 


0.380 


83 


.94 


52.34 


169. 


.2 





.407 


0.380 


-0.06 


83.89 


7.80 


137.46 


55.11 


158. 


.8 





.421 


0.400 


85 


.69 


54.82 


158. 


6 





,421 


0.401 


-0.05 


85.65 


8.30 


151.11 


57.12 


149 


.1 





.431 


0.416 


87 


.03 


56.85 


149, 





0. 


.431 


0.417 


-0.04 


86.99 


8.80 


165.13 


58.77 


140. 


.3 





.440 


0.429 


88 


.05 


58.52 


140, 


1 


0. 


,440 


0.429 


-0.02 


88.03 


9.30 


179.53 


60.15 


132 


.4 





.447 


0.439 


88 


.84 


59.89 


131. 


8 


0. 


.447 


0.439 


0.04 


88.88 
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